Genomic and Glycolytic Entropy Are Reliable Radiogenomic Heterogeneity Biomarkers for Non-Small Cell Lung Cancer.

Radiogenomic heterogeneity features in 18F-fluorodeoxyglucose positron emission tomography (18F-FDG PET) have become popular in non-small cell lung cancer (NSCLC) research. However, the reliabilities of genomic heterogeneity features and of PET-based glycolytic features in different image matrix sizes have yet to be thoroughly tested. We conducted a prospective study with 46 NSCLC patients to assess the intra-class correlation coefficient (ICC) of different genomic heterogeneity features. We also tested the ICC of PET-based heterogeneity features from different image matrix sizes. The association of radiogenomic features with clinical data was also examined. The entropy-based genomic heterogeneity feature (ICC = 0.736) is more reliable than the median-based feature (ICC = -0.416). The PET-based glycolytic entropy was insensitive to image matrix size change (ICC = 0.958) and remained reliable in tumors with a metabolic volume of <10 mL (ICC = 0.894). The glycolytic entropy is also significantly associated with advanced cancer stages (p = 0.011). We conclude that the entropy-based radiogenomic features are reliable and may serve as ideal biomarkers for research and further clinical use for NSCLC.

Elucidate multidimensionality of type 1 diabetes mellitus heterogeneity by multifaceted information.

Type 1 diabetes (T1D) is an autoimmune disease. Different factors, including genetics and viruses may contribute to T1D, but the causes of T1D are not fully known, and there is currently no cure. The advent of high-throughput technologies has revolutionized the field of medicine and biology, and analysis of multi-source data along with clinical information has brought a better understanding of the mechanisms behind disease pathogenesis. The aim of this work was the development of a data repository linking clinical information and interactome studies in T1D. To address this goal, we analyzed the electronic health records and online databases of genes, proteins, miRNAs, and pathways to have a global view of T1D. There were common comorbid diseases such as anemia, hypertension, vitreous diseases, renal diseases, and atherosclerosis in the phenotypic disease networks. In the protein-protein interaction network, CASP3 and TNF were date-hub proteins involved in several pathways. Moreover, CTNNB1, IGF1R, and STAT3 were hub proteins, whereas miR-155-5p, miR-34a-5p, miR-23-3p, and miR-20a-5p were hub miRNAs in the gene-miRNA interaction network. Multiple levels of information including genetic, protein, miRNA and clinical data resulted in multiple results, which suggests the complementarity of multiple sources. With the integration of multifaceted information, it will shed light on the mechanisms underlying T1D; the provided data and repository has utility in understanding phenotypic disease networks for the potential development of comorbidities in T1D patients as well as the clues for further research on T1D comorbidities.

Construction and Analysis of Protein-Protein Interaction Network of Heroin Use Disorder.

Heroin use disorder (HUD) is a complex disease resulting from interactions among genetic and other factors (e.g., environmental factors). The mechanism of HUD development remains unknown. Newly developed network medicine tools provide a platform for exploring complex diseases at the system level. This study proposes that protein-protein interactions (PPIs), particularly those among proteins encoded by casual or susceptibility genes, are extremely crucial for HUD development. The giant component of our constructed PPI network comprised 111 nodes with 553 edges, including 16 proteins with large degree (k) or high betweenness centrality (BC), which were further identified as the backbone of the network. JUN with the largest degree was suggested to be central to the PPI network associated with HUD. Moreover, PCK1 with the highest BC and MAPK14 with the secondary largest degree and 9th highest BC might be involved in the development HUD and other substance diseases.

Genetic signatures of heroin addiction.

Heroin addiction is a complex psychiatric disorder with a chronic course and a high relapse rate, which results from the interaction between genetic and environmental factors. Heroin addiction has a substantial heritability in its etiology; hence, identification of individuals with a high genetic propensity to heroin addiction may help prevent the occurrence and relapse of heroin addiction and its complications. The study aimed to identify a small set of genetic signatures that may reliably predict the individuals with a high genetic propensity to heroin addiction. We first measured the transcript level of 13 genes (RASA1, PRKCB, PDK1, JUN, CEBPG, CD74, CEBPB, AUTS2, ENO2, IMPDH2, HAT1, MBD1, and RGS3) in lymphoblastoid cell lines in a sample of 124 male heroin addicts and 124 male control subjects using real-time quantitative PCR. Seven genes (PRKCB, PDK1, JUN, CEBPG, CEBPB, ENO2, and HAT1) showed significant differential expression between the 2 groups. Further analysis using 3 statistical methods including logistic regression analysis, support vector machine learning analysis, and a computer software BIASLESS revealed that a set of 4 genes (JUN, CEBPB, PRKCB, ENO2, or CEBPG) could predict the diagnosis of heroin addiction with the accuracy rate around 85% in our dataset. Our findings support the idea that it is possible to identify genetic signatures of heroin addiction using a small set of expressed genes. However, the study can only be considered as a proof-of-concept study. As the establishment of lymphoblastoid cell line is a laborious and lengthy process, it would be more practical in clinical settings to identify genetic signatures for heroin addiction directly from peripheral blood cells in the future study.

Association between depressive disorders and risk of breast cancer recurrence after curative surgery.

The aim of the study was to investigate the association between depressive disorders and risk of tumor recurrence in patients with breast cancer after curative surgery.A nationwide cohort study between January 2001 and December 2007 was conducted. Data were taken from the Taiwan National Health Insurance Research Database. Among 30,659 newly diagnosed breast cancer patients, we identified 1147 breast cancer patients with depressive disorders and 2294 matched breast cancer patients without depressive disorders, both of whom received curative breast surgery between January 2003 and December 2007.The risk of first tumor recurrence was compared between patients who developed depressive disorders after breast surgery (depressive disorder cohort, n = 1147) and matched patients who did not develop depressive disorders (matched nondepressive disorder cohort, n = 2294). Cumulative incidences and hazard ratios (HRs) were calculated after adjusting for competing mortality.The depressive disorder cohort had a higher rate of recurrence when compared with the matched nondepressive disorder cohort (17.1% vs 12.5%; P < .001). The Kaplan-Meier analysis revealed a predisposition of patients with depressive disorders to suffer from recurrence (log-rank test, P < .001). After multivariate adjustment, the HR for subsequent recurrence among the depressive disorder cohort was 1.373 (95% confidence interval 1.098-1.716, P = 0.005). Moreover, the depressive disorder cohort had higher risk of overall mortality even though not significant after adjusted (adjusted HR 1.271, 95% confidence interval 0.930-1.737, P = 0.132).Depressive disorder was associated with a higher risk of breast cancer recurrence among patients after curative breast surgery.

Genotyping of human neutrophil antigens (HNA) from whole genome sequencing data.

BACKGROUND: Neutrophil antigens are involved in a variety of clinical conditions including transfusion-related acute lung injury (TRALI) and other transfusion-related diseases. Recently, there are five characterized groups of human neutrophil antigen (HNA) systems, the HNA1 to 5. Characterization of all neutrophil antigens from whole genome sequencing (WGS) data may be accomplished for revealing complete genotyping formats of neutrophil antigens collectively at genome level with molecular variations which may respectively be revealed with available genotyping techniques for neutrophil antigens conventionally. RESULTS: We developed a computing method for the genotyping of human neutrophil antigens. Six samples from two families, available from the 1000 Genomes projects, were used for a HNA typing test. There are 500 ~ 3000 reads per sample filtered from the adopted human WGS datasets in order for identifying single nucleotide polymorphisms (SNPs) of neutrophil antigens. The visualization of read alignment shows that the yield reads from WGS dataset are enough to cover all of the SNP loci for the antigen system: HNA1, HNA3, HNA4 and HNA5. Consequently, our implemented Bioinformatics tool successfully revealed HNA types on all of the six samples including sequence-based typing (SBT) as well as PCR sequence-specific oligonucleotide probes (SSOP), PCR sequence-specific primers (SSP) and PCR restriction fragment length polymorphism (RFLP) along with parentage possibility. CONCLUSIONS: The next-generation sequencing technology strives to deliver affordable and non-biased sequencing results, hence the complete genotyping formats of HNA may be reported collectively from mining the output data of WGS. The study shows the feasibility of HNA genotyping through new WGS technologies. Our proposed algorithmic methodology is implemented in a HNATyping software package with user's guide available to the public at http://sourceforge.net/projects/hnatyping/.

Development of a questionnaire on everyday navigational ability to assess topographical disorientation in Alzheimer's disease.

We developed a Questionnaire on Everyday Navigational Ability (QuENA) to detect topographical disorientation (TD) in patients with Alzheimer's disease (PwAD). In the QuENA, 3 items were designed to assess landmark agnosia, 2 for egocentric disorientation, 3 for heading disorientation, and 2 for inattention. The PwAD and their caregivers rated QuENA according to which TD symptoms would occur. Regarding the construct validity, confirmatory factor analysis showed that the caregiver version of the QuENA fits the proposed TD model well but the patient version does not. Regarding the internal consistency, the Cronbach's α for the caregiver version was 0.91 and that for the patient version was 0.87. A discrepancy existed between the appraisal of navigational abilities by PwAD and by caregivers, and it was correlated with the number of getting lost (GL) events. The caregiver version of QuENA is a feasible, reliable, and valid instrument to assess TD and it also discriminates well between the PwAD with GL and those without.

Associations between HIV and human pathways revealed by protein-protein interactions and correlated gene expression profiles.

BACKGROUND: AIDS is one of the most devastating diseases in human history. Decades of studies have revealed host factors required for HIV infection, indicating that HIV exploits host processes for its own purposes. HIV infection leads to AIDS as well as various comorbidities. The associations between HIV and human pathways and diseases may reveal non-obvious relationships between HIV and non-HIV-defining diseases. PRINCIPAL FINDINGS: Human biological pathways were evaluated and statistically compared against the presence of HIV host factor related genes. All of the obtained scores comparing HIV targeted genes and biological pathways were ranked. Different rank results based on overlapping genes, recovered virus-host interactions, co-expressed genes, and common interactions in human protein-protein interaction networks were obtained. Correlations between rankings suggested that these measures yielded diverse rankings. Rank combination of these ranks led to a final ranking of HIV-associated pathways, which revealed that HIV is associated with immune cell-related pathways and several cancer-related pathways. The proposed method is also applicable to the evaluation of associations between other pathogens and human pathways and diseases. CONCLUSIONS: Our results suggest that HIV infection shares common molecular mechanisms with certain signaling pathways and cancers. Interference in apoptosis pathways and the long-term suppression of immune system functions by HIV infection might contribute to tumorigenesis. Relationships between HIV infection and human pathways of disease may aid in the identification of common drug targets for viral infections and other diseases.

A quantitative analysis of monochromaticity in genetic interaction networks.

BACKGROUND: A genetic interaction refers to the deviation of phenotypes from the expected when perturbing two genes simultaneously. Studying genetic interactions help clarify relationships between genes, such as compensation and masking, and identify gene groups of functional modules. Recently, several genome-scale experiments for measuring quantitative (positive and negative) genetic interactions have been conducted. The results revealed that genes in the same module usually interact with each other in a consistent way (pure positive or negative); this phenomenon was designated as monochromaticity. Monochromaticity might be the underlying principle that can be utilized to unveil the modularity of cellular networks. However, no appropriate quantitative measurement for this phenomenon has been proposed. RESULTS: In this study, we propose the monochromatic index (MCI), which is able to quantitatively evaluate the monochromaticity of potential functional modules of genes, and the MCI was used to study genetic landscapes in different cellular subsystems. We demonstrated that MCI not only amend the deficiencies of MP-score but also properly incorporate the background effect. The results showed that not only within-complex but also between-complex connections present significant monochromatic tendency. Furthermore, we also found that significantly higher proportion of protein complexes are connected by negative genetic interactions in metabolic network, while transcription and translation system adopts relatively even number of positive and negative genetic interactions to link protein complexes. CONCLUSION: In summary, we demonstrate that MCI improves deficiencies suffered by MP-score, and can be used to evaluate monochromaticity in a quantitative manner. In addition, it also helps to unveil features of genetic landscapes in different cellular subsystems. Moreover, MCI can be easily applied to data produced by different types of genetic interaction methodologies such as Synthetic Genetic Array (SGA), and epistatic miniarray profile (E-MAP).

Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme.

BACKGROUND: The development of microarrays permits us to monitor transcriptomes on a genome-wide scale. To validate microarray measurements, quantitative-real time-reverse transcription PCR (Q-RT-PCR) is one of the most robust and commonly used approaches. The new challenge in gene quantification analysis is how to explicitly incorporate statistical estimation in such studies. In the realm of statistical analysis, the various available methods of the probe level normalization for microarray analysis may result in distinctly different target selections and variation in the scores for the correlation between microarray and Q-RT-PCR. Moreover, it remains a major challenge to identify a proper internal control for Q-RT-PCR when confirming microarray measurements. RESULTS: Sixty-six Affymetrix microarray slides using lung adenocarcinoma tissue RNAs were analyzed by a statistical re-sampling method in order to detect genes with minimal variation in gene expression. By this approach, we identified DDX5 as a novel internal control for Q-RT-PCR. Twenty-three genes, which were differentially expressed between adjacent normal and tumor samples, were selected and analyzed using 24 paired lung adenocarcinoma samples by Q-RT-PCR using two internal controls, DDX5 and GAPDH. The percentage correlation between Q-RT-PCR and microarray were 70% and 48% by using DDX5 and GAPDH as internal controls, respectively. CONCLUSION: Together, these quantification strategies for Q-RT-PCR data processing procedure, which focused on minimal variation, ought to significantly facilitate internal control evaluation and selection for Q-RT-PCR when corroborating microarray data.

Effect of Epstein-Barr virus infection on global gene expression in nasopharyngeal carcinoma.

It was proposed that Epstein-Barr virus (EBV) is closely associated with nasopharyngeal carcinoma (NPC); however, the molecular mechanisms involved in the effect of EBV on NPC host genes have not yet been well defined. For this study, two sets of microarray experiments, NPC (EBV-free) vs normal epithelial cells and EBV(+) vs EBV(-) NPC arrays, were analyzed and the datasets were cross-compared to identify any correlation between gene clusters involved in EBV targeting and the NPC host gene expression profiles. Statistical analysis revealed that EBV seems to have a preference for targeting more genes from the differentially expressed group in NPC cells than those from the ubiquitously expressed group. Furthermore, this trend is also reflected in log ratios where the EBV target genes of the differentially expressed group origin showed greater log ratios than genes with an origin from the ubiquitously expressed NPC group. Taken together, the genome-wide comparative scanning of EBV and NPC transcriptomes has successfully demonstrated that EBV infection has an intensifying effect on the signals involved in NPC gene expression both in breadth (the majority of the genes) and in depth (greater log ratios).

A stochastic differential equation model for quantifying transcriptional regulatory network in Saccharomyces cerevisiae.

MOTIVATION: The explosion of microarray studies has promised to shed light on the temporal expression patterns of thousands of genes simultaneously. However, available methods are far from adequate in efficiently extracting useful information to aid in a greater understanding of transcriptional regulatory network. Biological systems have been modeled as dynamic systems for a long history, such as genetic networks and cell regulatory network. This study evaluated if the stochastic differential equation (SDE), which is prominent for modeling dynamic diffusion process originating from the irregular Brownian motion, can be applied in modeling the transcriptional regulatory network in Saccharomyces cerevisiae. RESULTS: To model the time-continuous gene-expression datasets, a model of SDE is applied to depict irregular patterns. Our goal is to fit a generalized linear model by combining putative regulators to estimate the transcriptional pattern of a target gene. Goodness-of-fit is evaluated by log-likelihood and Akaike Information Criterion. Moreover, estimations of the contribution of regulators and inference of transcriptional pattern are implemented by statistical approaches. Our SDE model is basic but the test results agree well with the observed dynamic expression patterns. It implies that advanced SDE model might be perfectly suited to portray transcriptional regulatory networks. AVAILABILITY: The R code is available on request. CONTACT: cykao@csie.ntu.edu.tw SUPPLEMENTARY INFORMATION: http://www.csie.ntu.edu.tw/~b89x035/yeast/